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I. DESCRIPTION OF THE DIFPUSLON MODEL 


A. INTRODUCTION 


The main objective of this thesis is to fini and test an 
experimental mathematical model which predicts the prop- 
Moi ity of detecting a two-dimensional tacget by “a 
Stationary searcher. This model will be shown to provide an 
upper bound for the probability of detection by a stationary 


Beatecher OL a target Conducting a "random tour” ( Ref. 1] 


Bee DESCRIPTION OF DIFFUSLON MODEL 


lie We SSee Seige Mieresiesloyel 
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The searcher is assumed to be located in the center 
of a sguare search region of area A. This location is held 
fixed during the search period. The searcher has a detec- 
mren Capability over a disk of radius &£&. The detection 
probabiiity of a target inside of this disk is 1 and outside 
is 0. The searcher thus has a "cookie— cutter" sensor with 


detection range FR. [Ref. 2} 


The EdegeulS (St art in g poSition is uniforaly 


distributed over the sguare search region A. 


3- Motion of the Target 


In our diffusion model, the target moves randonly 
over the area A as a diffusing pacticle which ceflects off 
tre area boundaries. The diffusion constant is D, which has 
dimensions of area per unit of time. In any time interval 


of length At that does not contain a boundary reflection, 


components of the target's position on the X and Y axes 
suffer increments which are independent o£ all previous 
increments and which are each distributed normally wath mean 
QO and variance DAt. 

Still ignoring boundary ~er2ceuc, the diritusies 
assumption results in the target's location at time t having 
a circular bivariate normal probability distribution wakes 
mean of the starting position and variance of Dt. 

Thus the probability density of the tacget's loeae 


tVOn dt  ermesarrs 


i 2 
| (X-U,) + (y-Uy) 
£ =e oe 
E {ep ye =) 27T Dt exp 2D0t (7.1) 
where (Uy .U,) dso) thestar Gets rsie a remem oe oem ene re Adding 


the effects of boundary reflection significantly complicates 
tne -Caleculat1 on Oummrs (> and Leads to the necessity of 


uSing Simulation to attack this problen. 
4. Detection 


Detection occurs whenever the target enters the 


searcher circular detection disk which has a radius R. 


C. DIFFUSION SIMULATION MODEL (DIFSIM) 


A Monte Carlo Simulation computer model (DIFSIM) is used 
to genarate detection probabilities for this diffusion 
model. This projram is written in FORTRAN and designed for 
use at the Naval Postgraduate School(NPGS). It uses the new 
version of the NPSS Random Number Genarator Package, called 
LLRANDOMII in order to generate Jniform and Normal randon 
numbers. 
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1. Inputs 


—_— = 


e Area size, A, in scuare nautical niles. 


e Dirfusion constant, D, in sguare nautical miles per 


mouse 
e Padius of detection jisk, R, in nautical miles. 
e Numter of replications. 
® Detection period as an Eour. 
e Time increment, At, for each discrete step in minutes. 


2. Functioning of Progran 


The initial target position is selested irom a 
bivariate uniform distribution over the search region A. 
Subsequent target positions are determined by a discrete 
approximation of the diffusion. tie make the following defi- 


nitions, 


X=x component of current location 

Y=y component of current Location 

A*‘=x component of cew Location at the end 
of time increment At 

Y*'=y component of new Location at the end 


of time increment At 
Then, 
k1=x+ © (pat)? 
eer a eae 


where 6 and OQ) are drawn independently from a standard normal 
ssc ribution. 
In this model a 5 minute At is used. Different time 


increments, varying from 1 minute up to 15 minutes, have 


11 


been tested and 5 minutes has been accepted as 4a good value. 
For Smaller time increments the Simulation projcram too mec 
lony in computer execution tine. As Shown in FPiguiurce wee 
thece is ho significant dirference in probability curves 


between 5 minute increments and smaller time increments. 
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When the target encounters a boundary, i rerlectaies 
is made to keep the target insiie the search area. Ine 


target's Y¥' position after a reflection 1S giver as toute 


Yt<20 => Y' becomes --Y'! 


Y'> a => Y' becomes 2Za-Y' 


where a is the length of a side of tne Square Search area A. 

The target reflects in the X direction in a Similar tanner. 
Detection occurs whenever the target enters the 

detection disk. This event can be defined analytically as 


TOL. Ows: 


a a 
SoS See) SS Be 


For each time t, the Simulation output is the ratio 

No/N, where 
Njo=humber of replications jiviny a detection by tine t, 
N,=total number of tionte Carlo Replications used in the 


Simulation. 
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II. CCNRECTION BETWEEN RANDOH TOUR AND DIFFUSLON MODEL 


A. DESCRIPTION OF RANDOM 4 CURS MODEL 


In the random tour model considered here, the target is 
assumed to move at a constant speed and td make course 
changes at randon times. Each new course iS drawn crom 
Ulitf OLN Gistr@outLon on, Oy iin. Ihe lengths »9f tne time 
between course CaangeS iS exponentially distributed with 
mean t/A. 

An analytic expression for the probability density of 
the target's positicn after a random tour of length | tive 
derived in {Ref. 1] Given the taryjet's initial position at 
the origin of a two-dimensional coordinate systen, this 


expressicn is 


g(x,y,t) STWR jee (-At ( 1- /7-$)) rexp™ (S(¢-1) (2. 1) 


V=target speed (nautical mile in per hkour) 
A=course Change rate (nrs ) 


t=time (hrs) 


a a 
g X +yYy 
(vt) 
The Dirac §_function component orf g9(x,y,t) arises ese 


Ene SetaGteeuac. with probability e the targst does nor 
finish the first step by time t. In other words, the targes 
Makes no course change through time t. Therefore, its prob- 
ability mass is concentrated on the boundary of a disk of 


radius vt and centered at the origin. 


iG 


PCC eee s ia brn PhO TOUR Alo DirFUSIOH MODELS 
eo Geena tse «6 OL UTR TAL, kt Lise. th: 
Percoll at wime tL. I nen 


Re=xe+ Ye 


Sue expression for E£[{R2] cor the random touc model was 


meer ved In [Refs 3] as collouvs: 


siaeye | f cxeey? g(*,y,t) dxdy eas) 


O<x2+y2< (vt) 2 


BeoStituting (2. 1) Sey eoeen (en iw ted SrOnrnn hy “LO polar 


coordinates, we obtain 


ar WA 





§j—=-  - 1] + aa rt /1-{— -drdé 
te Le " Vt 
: | 
Seating x=r/vt, we have 
z 


E[R?] = eet (Vt)2 fx? ¢ 6(x-1) 


0 
ae oes 
RT 


= “vey? 1+ ops exp 4 At fr dx (2.3) 


tox ( ne [Te] dx 


NO 





le 


. : . . = - ee = ey = 
To “1 riopw the Ln vegieeO ni eae ie aie t d=) ia ; oe! 


CGE ey: aes 


1 
3 
tr | ——— eXxD dt} l-x? dx 
1-x? 
o { 


= et (ete) exp ( rev) du 


0 At 
wae =a 7 y e? dy 
Cre 
0 
ht 
_ _2e Cel = ea eee (2.4) 


From (222) sand (2.4) Sweetie ne hawe 


2 _ 
E[R?] = 2 ee Ore sie ile ate Seay 
pee 
_ tye f,. tee 
» A Sag) 


In the calculation of Ej, R2] for diffusion nodetiieaa 
difiusion constant D, the tarcet's initial positionmmee 
assumed to te at the origin of a two-dimensional cooramiae. 


system. We have 


Ei R22 j-h Xe +e |= sexe yeas (2.6) 


es 


since X and Y are independent and uncorrelated, they have 


Mopial Darginal distributions. So, 
X~ i (0,Dt) 


Y~N(0,Dt) 


an Gd. 
gf X@j=Vac[ X J+ (2[X]) 2=Dt+0=Dt 
E{ Y2 J=Vacl Y J+ (E[¥)}) 2=Dt+0=Dt 


If we sukstitute these Ej X2?}] and E, Ye] values into (2.6), we 


get 
E( Re J=EL Xe J +EL Ye J=Dt+Dt=2dt 2a) 


As described in {Rez. 3] tor tne random tour model, ast 
goes to infinity the Central Limit Theorem rejuires that 
Mee y,e) becomes asymptotically circular bivariate norma. 
With mean Pp =9 and variance g?=Vt/A. This result can be 
Serained by tising the formula (2.5) and letting t go to 
mea nity « 

"av't 


| as a cas (2.5) 


aV't |~e 
At A 


Lim E{ R@j=Lin y ie 


t--> co t--> co 


Ey wcoOllbdriiiecneseqiatrons {72.7}. and (2.3), it 1S Seen 
that as t becomes large, a random tour can be approxianated 


aed GlLTuUSi0n with diffusion constant 
=a (23) 


To examine the relationship between a random tour and 
its “eguivalent" diffusion, two simulations were used. 


Example results of these two Simulation programs, DIFSIM anc 


i 


PASS are displayed in Figure 2.1. Both progjgrags are Monte 
Cario sé€arch sSimtlattons. For PASS (Passive AcousStume 
Submarine Simulation), [Ref. 4] the target motion model isa 
random tour. In DIFSIM (Diffusion Simuiation), the tapeeg 
moves as a diffusing particle. In botn cases, the sSearcner 
is stationary at the center of a 100 rm x 100 hn Seater 
area and has a cookie cutter sensor with detection range 
15nm. The target reflects off the aseavounrdia cies. The 
initial target position was selected uniformly over the 
search area, and each replication or the Simulation was 
ailowed to continue until the target moved within distance 
15 nm of the searcher. 

To generate the results shown in Figure 2.1, the 
foliowing 5 different pairs of A and VY are used 45 a raceme 
course change and speed of target in tke PASS simulation 


model. 


v 5.0\ /10.0 15.0 20.0 25.0 1 O/ AE 
A 0.25), \tl0/ © eee aay en eres 1/hr 


If we use equation (2.9), we na y get an "eguivalent" difftu- 
sion constant 100 nm@/hr, by using each ditferent  (A,V) 
ba eee Thus 100 nm2/hr is used as a diffusion constant in 
the DIFSIM in order to get an “equivalent"™ ditfiusion modeum 

As demonstrated, detection probabilities are asymptoti- 
cally very close to each other as t increases. But during 
the early search hours, datection srobabilities for a dizfu= 
Sion model are higher than the probabilities which are 
generated by the random tour model. 

If we recall the the equations (2.5) and (2.8) we wae 
see that 


2 
E[ R2 ]S---<---=20Dt 
A 
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Thus the approximation of E[R2] for the candom tour 
model by using the diffusion model always leads to an 
OVERESTIMATE of EL R2). 

In the diizfusion model we may expect that the taryet 
Will move a greater distance from the origin taan does the 
Saeget in the “equivalent” random tour model. It is there- 
fore reasonable to expect that the diffusing target wili 
will <a 


Corn 


se 


a 
hi 


ture has been supported by further Simulation testing; also 


ap 


Gt 


encounter a stationary searcher aore guickiy 


ano” 
Sie 


bs 
UV) 
Ld. 


fanget Conducting an "equivalent" random tour. 


Supported is the fact that the two processes are asymptoti- 
Gaity identical. 
AN experinental analytical model will be constructed for 


the diffusion model in the next section. 
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III. ESTABLISHMENT OF THE ANALYTICAL MODEL 


Poe tiisechap ter, dnd wWhenesimuiiti On results Erom DIFSIM, 
aimee xpeLbimental analytical model will be constructed for 
predicting the probability of a diffusion target entering a 


stationary disk by scme time t. 


A. ASSUMPTIONS 


The followiny assumptions are used in our nodel 
eTne searcher is fixed at the center of a sjuare search 


region of area A. 


eThe searcher has a detection capability overc a disk of 
BAd User with da Getectaon probability of 1 within the 


dask and 0 outside. (ise, n a  COOKre-CULteL A2etector). 


Sine tabget'sS Starting position is uniformly distributed 


over a search region A. 


eThe target moves randomly over the area Aas a diffusion 


naeereleuwitieciciiicton constant D- 
eThe target reflects off the area boundaries. 


eA target can be detected only once by the searcher 


Bee CLASSIFICATION OF THE YARIABLES 
The variables in our model are, 
Srirerpaaereece 2 On proba bility, P. 
eSearch area, A (nm?). 
eSearcher detection disk radius, R (nm). 


eTarget diffusicn constant, D (nm?/hr). 


a3 


eTine, t« (ha). 
The first detection probability, P, is the depsndene yaa 
able. The renaining variables A,k,D and t ars inaependent 


Ea se 


p= £(A,R,D,t) 


C. ‘CONSTRUCTION OF THES Ope 


Figure 3.1 shows four piots of the) propactiity sou 


Oo ® 


target detection Ly time t as estimated ty the Honte Carl 
SIMULATION PLE Sis If we look at these curves, we wiil 
ofserve that all of them have an increasing trend and tie 
acproach 1 asymptotically. Pt alsoea9e3 as eae ee 
derivative must be negative everywhere. Figure 3220 pileu sae 
Minus the same data on a log scale. The fact that these 
plots are very neariy linear suygests the following func- 


tiojia Peroni se Tt) 
errno ca as (3. 1) 


where a and gp are determined by the problem parameters &,A 
aid aoe After conducting 23 sevarate simulations ¥aed 
different values of R,A and D, the author is convinced that 
the form of P(t) 1S approximately exponential. This thesis 
attempts to establish values of a and B87 ds functtomeo mere 


and D, to ailow eguation (3.1) to be a reasonable estimate. 
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In our model, we know that the target starting posi- 
tion is uniformly distributed over the search area a, and 
the searcher has a detection capability over a disk wuich 
covers an area [fTP2 with probabilty 1. Soy w2 May expect 
that at the beginning of the Search, i.e., when t=0, detec- 
tion probability will be egual to/TrRe/A. 


If we subsititute t=0 in eyuation (3.1), we get 


P(O)= 1-& 
Then 


E(t=0) een x 2) 


which implies that 
a = 1--—————— sien 3) 


yo ou omede. for B 


ie ooodeimidi imine hide ait andepenisit variables 
R,D,A and t, and will therefore be more complex. “@ Will 
study each independent variable Separately in order to 
Simplify the problem (We wili change one variable while 
holding the others fixed.) Tne relationship between B an 4 
these variables will be estimated. Eventuaily we will 


combine these submodels for a final submodel. 


a. The Relationship Between Diffusing Constant PD 
and B 


For this case, area size A and craljius & were 
held fixed at 10000 nm2 10 nao. respectively. piriuston 
constant D was varied between 20 nm?@/nr and 30090 mm?/nr The 


Simulation results are displayed in Figures 3.3. 


ar 


By using the least syuares estimation method on 
Inf 1-P(t)] data values, a best fit 6 was obtained Tore ee 


diffusioa constant D. 
These B values are plotted against the  jcoree. 


sponding D values in Figure 3.4. They fall approximately on 
destrdal diese ioe. Again, ody using the least syuare estita- 


tion method, the following linear equation was outained 


B =0.003038+0.000205 D (3.4) 
™~0.000205 b 
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bk. The Relationship Between Area Size A ana B 


This time the diffusion constant D aid detectaaa 
radius R were held at 100 nm?/hr and 10 sm c2spectively. 
The area size A was varied between 1000nm* and 20000 nme. 
The Simulation results and the least square estimation 
results for the log transformed data values are displayed in 
Figure 3.5 and and Table II 

Figure 3.6 shows a plot of the best fit p 
against area size A. These points fit very closely with the 


power function 
B =O0.77A71-49 


(A least squares procedure on the log transforned data was 
used to determine the values 0.77 and -1.49). To achieve 
moce natural final units, the model is modified slightly as 
follows: 


B =0.77A-1-5 (B25) 
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c. The Relationship Between Detection Radius R and 


B 


Again, the same procedures were applied for this 
case. R was varied between 1 and 30 nm, while D and A were 
held fixed at 100 nm*/hr and10000 nm? respectively. 

The simulation results and least sgiare estimna- 
tion results are displayed in Figure 3.7 and Table Ii 13a 
Figure 3.8, the scatter plot of estimated f> values shows a 
linear relationship) betwee te hanes. 


Least square eStimation for this lin2 is 


 =-0.00105+0. 00278R (3.6) 
Bx 0.00278R 
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We can summarize tnese obServations is rollows-. 
po D 


ee A 715 


Bo k 


where "«%" means ‘is proportional to" waich sugg2sts that 


RD 
Pa 


(3-7) 
for the rroper value of K. To estimate K is the finai rodel 


building step. 
dad. -Fstimaticn of the Cociereteontce.s 


We can calculate the K value for each SinulatiGcz 


with the expression 


rs. 
PA (3:3) 
RD 

where B is the "best fit" value for that similation ie 
Then with these sample &£ values we may find best overall 
estimate. 

In addition t9 56 simulations already conpletcay 
25 additional simlations were conducted to produce a total 
of 81 sample kK values. The histogram and the statistical 


table vaiues for this data are displayed an Pigarce 322e 
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Tf we recall the ¢€ juation {351)7 (3.2) fae eee 


Pete) = ] - ue 
ao = ] - is 
A 
= Se 
Ais 


Substituting these a and @ values in eqWaeaaem 
(25) we derived our final analytical model for first 


detection probabilities as follows 


RD 
aR? = Ke sa: eee 
2 = -_ —— .9 
Por) it Gal a Vere fers (3.9) 
where P(t) refers to the probability of Eirst deteeumem 


OCCUrS ONvOmEDeLOLe stizeat 


Dey VERIEICAIION Of Sint ee DEE 


Ee = SESS —_ = SS — ae eS 


From equation (3.9) we see that 


ie ae 
a 


must be dimensiouless. This inplies that the coefficienmeme 
must be dimensionless. (If we had set the power of A to 
1.49 versus 1.5, then the Jinension of K would be nm%2, not 


a td Cuba luni...) 
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Tf we hoid faxed all indevendent variables &,D,t and 
change the area size A, we will observe that as we increase 
the A, the probability 9f detection decreases and vice 
VerSa. This result iS demonstrated in Figure 3.10. AS we 
increase the size of tke search aréa, more area wiil be 
available for the target to eScape fron the Searcher. 
Therefore we may expect lower detection probabilities for 
larger search area sizes. 

A Similar SenSitivity analysis is applied tor the 
other indepenient variables &,D and t. If we look at the 
results displayed in Figure 3.10, we may observ2 that as we 
increase these variables, the detection prokabiiities 
increase Simultaneously. These results seem freasonabkie, 
because aS we increase the search time or detection radius, 
we may have more chanses to detect the target. <Aiso an 
increase in the D value neans that the targe Wait raves 
more distance during any time interval and will thus be nore 


likely to enter the detection disk. 
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3. Final Verification 


There was no actual data available fron real life 
observations. Therefore, the output of DIFSIM is used for 
final verification of the model. For this purpose, conbina- 
tions of the following independent variables wece used both 


miour dbdlytical model and aS input to DIFSi«. 
D= 40, 80, 140, 200 Nu*2/Hr 
Rec. 2) gee) Ogee On Zo Nm 


A= 4000, 8000, 12000, 16000 Nn*2 

The outputs are displayed in Figure 3.,11, Bi lo 
Bad 3. It is observed that the simulation and model prob- 
abilities are generally very close to each other. Only 
during early search hours do the simulation curves sometimes 
go akove the model curve. This means that our model 
predicts fewer detection than the diffusion simulation nodel 
during the early search hours. Since this difference is at 
Maximum .03 or .04, we conclude that the model provides a 
good fit, which gets better for larger time t. For smali 
t,the medel appears to underestimate the prosgability of 
detection. 
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ie ete eee Gc the 11k bait CUS et 
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i ee | Ore eG Ou ein 


Beeeetctethat the @irSih 7;eebabrlity of Uctection iy tipe t 
generally exceeds that predicted by the model suggests that 


tne instantaneous probability of detection produced by 


DIFSIM exceeds (Seo ledste sine Carly hours of the Siryla- 
Ory ) . To test this nypothesis, DIFSIN was cun with the 
Baramneters 


D=100 om@ /hr 
R=10 Do 
A=10000 nm 


[eeu «Of Cach 5 hour period between 0 and 80 hours, tke 


jeast squares best fit was obtained. (That 1S, the best 
tt 1 NG 
2 
ie lo ae ee 
A 


was Calculated). These Waites are PVilGtitcayin Figure 3. 14, 


and appear to approach from above the model value of 


mee exe LOO x10 E0267 
10000 


Pico eID thescdse at least, it appears that the irstanta- 
Meous probability of detection starts at some high value and 


decreases asyoptotically toa steady state value given by 


ug 


the medel presented here. This observation remains to be 


proved in general. 
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Ee CULULATIVE DETECTIONS 2 OSA cilie as eae. 


Since the following properties fOr general culmulaeae 
,roLability function hold £Or our tatheretical (ocala 
tion (3.9%), we tay assume that this model also represenocms 
Cumulative density Luneumeom cee) for first detection Ga 
flere 


These properties are 


Vee Liter fe) — Sa 








t—> & 
2 RO 
-K t 
Lin F(t)= Hs} 1-(1- dis ) Als 
A 
_ 4. fe eR ee 
= 4 ne 


2eF {t) iS a nondecredsing Lumet uon 


If we take the the first derivate of F(t) with 


respect to t, we get 


RD 
dPCE) 2) RD eee ce Se 
dt ie A 


Since this equation always has nornecative values, we 


fay assume that F(t) 1S @@nondeemeascinguenve: one 
3-F (t) aS a COnNtinuouSs = =i newer, 


Therefore we may def iheeedi mace come nag. 


ae 


eo eee ee fm et) 


F(lty = Je ree e Ale? Gar. i) 


Die Peete, Clim Chen ll TY SD eisiey FruNrcizol)L 


Meme mock Vetted tose ony wreOnotbuli ce densi ty LuRction 
ee), L(t), PDyeecti RG Ue te ween rSt derivati veo Edi with 
mes,ect to tite t. 


f(t) _ 


ue et 








eo) 





it ke a Coos (Ue eronwO toro we hay e¢ 


pee dt = (x82) (1-58) eat ay 


i TR ~K4AS 
A eC 
0 


SiS 


Since this integration 1S not egual to 7, the eGsiatvon ite 


doesn't represent a Proper density LWietion tone. 


C. EXPECTED FERST DETECTIONS ait 


Let F(t) denote the cumulative nondetection probalb i sige 
PInet Lou (ena: 


F(t)= ©(No detection ip to time aa) 
=1-P(T<gt) =1-F (t) 
RD 
= TUR” \ eee 
F (ts ( A }e" o (4.3) 





The following formula can be used to find the expected 


detection time §/ 7). 
oo 


zi [Fw dt (4.4) 


0 
if we substitute (4.3) in - 4) 


cms (| otis ne dt 


0 


ae Q2 RO Be —t 
i I A @ dt 


“a0 (| A 
Als TR “Ket ~ 
=KRO 2 “€ 
a aa 0 
~ KRO (I- A 
AN _AVITR 
ELT =o RD (4.5) 
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eee @ toe gee este detection times are displayed for 
Opt rerent difrusion constantS (9), area sizes (A) and detec- 


tion radiuses (R) in Figure 4.1. 


wen CONDITIONAL DETECTION PROBAEITLITY FUNCTIONS 


If we assume that there will be no detection at the 
beginning of the search period, we may derive tne following 


@ampaitionali cdf. 


Fo (t)=P (Detection up to time t/no det. at time Q) 


eer tO sth 
PLT 20] 


plo<Tét] 


Pao) (4.6) 


If we substitute t=0 in eguation (4.3), we get 


2 
ro (WS 10) = IPO) (a-& (4.7) 
and 
P(O<T<t)= F(t) -P (t=0) 
2 
- fin (BE) okart|- OF (4.8) 


Bemuscing (4.7) and (4.8) as a dominator and numecatcr in the 


ecuation (4.6), we pave 
Fi(t)=- | oe pa Als t 
(4-9) 


5D 


This function (4.S) teapresents a cumulative probabalie, 


function of an exponential diStribatior with pdazamemem 


Ke By uSing this fact, csonditional expected Ea 
detection time E[I, ) ecammbe sdefiniea mas: ol tens: 
pio 
ELT, Jgsl L/No detec Blom sine aaa (4.10) 


Also, we can write conditional detection probability denaaaag 


cameron Eg (t) San tiesto trow nga tora 


RO 
RO a) cacretr t 
ft)2K a ae (4.11) 
If we compare equations (4.5) and (4.10), we will 


observe that, 


Ve) is LS 
py A ITR é A 
K RD KD SRD 
So fi) Tja56 aa This ineguality means that thea conditional 


first detection time is greater than the tuncondit lonal sea 
detection time. We can get this conclusion int vuitiveraaa 
thinking that we have an opportunity to detect the target at 
the beyinning of the search period in the unconditionai 
case. Thererore, for the unconditional case, w2 may expect 
an earlier first detection time than would be possible for 


Ena eS CconditiOmas cases 
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Figure 4.2 Corparison of F(t) and F,(t)- 
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ed 


Ve. THE RESULIS AND APPLICATION 


Be EXPECTED DETECTICN TIME FOR RANDOM TOUR MODEL 


As we showed in Chapter II, a diffusion target gives an 
upper bound detection probability for the "“eguivalent" 
random tour model. Therefore, the diffusion model expected 
detection time, which is estimated by eguation (4.5) should 
be a lower bound for random tour expected detection time. 


mat is, for 
Fi Tp j]=Diffusion Model Expected Detection Time, and 


FECT, J=Random Tour Model Expected Detection Time, we 


have 
Pe 7p 
tea — KO #7] 


Be ONE-DIMENSIONAL DIFFUSION MODEL 


In this thesis, two-dimensional diffusion motion was the 
rkasis for our model. The exponential type curves were used 
to estimate this model's outputs. Pree ad dist ON gO ume es 
model, the one-dimensional diffusion model is simulated by 
computer program DIFSIMN1. In this model, the target moves on 
a line segment L, according to diffusion constant D. The 
target's starting position is selected uniformly over this 
line ségment L. Detection occurs whenever the target hits 
the designated end point. The target reflects off the cther 
end pcint of line segment. 

The results of different Simulation results are 
displayed in Figure 5.1. Exponential curves, which were 


obtained by using the least squares estimation method, were 


a3 


used to estimate the outputs of one-dimensional dirtfusion 
model as we did for two-dimensional case. 

So we may expect that, for three or more diwensional 
diffusion models, we may use exponential type curves which 
are generated by a different set of parameters. For three- 
dimensional case, these parameters can be defined as 


follows: 
V=Volume of the cubical search space. 
R=Radius of the cookie-cutter detection Sphere. 
D=Difiusion constant. 


T=Detection time. 


C. APPLICATIONS OF THE MODEL 


Our model can also be used to estimate the final detec- 
tion prokability of a system which includes more than one 
independent sensor. As an exampie, we may use the following 
scenario: 

we want to use n sonobouysS in order t> aetect a 
diffusing target in an area A. Each sonobouy has a 
cookie_cutter detection capaktility over a disk with radius §; 
; Each sonobouy will be independetly located on the center 
of a square subsearch area A-and operatedwe for a tame pence 
i If we make the following assumptions, we may use egua- 
tion (3.9) to estimate the overall detection probability of 
this sonobouy pattern at the end of the search period t. 


ial 


where Ris the effective detection radius. 
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tila wo SOndiuded SO rg1iveswUsS an ~upser bound detection 
probability for a random tour model which moves in tne sane 


systen. 
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APPENDIX A 
DIFSIM COMPUTER PROGRAM 


in Order to give access to tae logic used in bduildias 
the Simulation models DIFSIM and DIFSI41, a complete progran 
Meseing iS imcluded anh this Appendix fEollowing the list of 


variaktles used in the Simulation nodels. 


UST OF VARLABLES 
REP =Number of replications. 
MAX =Detection period as a minute. 
k =Radius of detection disk in nautical niles. 


DEF =Datfusion constant Din Square nautical miles per 


MGW 


SIDE=The length of the scuare Search area side in 


nautical miies. 
AREA=Area size A in syuare nauticai niles. 


INC =Time increment t for each discrete st2p in 


minutes. 
FROB=Probability of detection 
POSX=X component of target position. 
POSY=Y component of target position. 


DIST=The distance between the target location and the 


CentenswoLr the detection disk: 
TX1 =Seed number for uniform random number. 
imo) —=Secdeeumber £LOr Standard normal random variabie. 
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